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We study a 3/4-filled two-dimensional extended Hubbard model under the fluctuation-exchange approx- 
imation. We find that this model undergoes phase separation in a region of nonzero temperature, where 
the quantum critical phenomenon of charge ordering dominates. By considering the long-range Coulomb 
interaction that frustrates this phase separation, we present a mechanism for generating a glassy state in a 
quasi-two-dimensional organic conductor on the verge of charge ordering. 

KEYWORDS: organic conductors, charge ordering, phase separation, glassiness, quantum critical phenom- 
ena 



Organic conductors are providing a fruitful stage for study- 

[ ing exchange-correlation effects between electrons on a crys- 
tal lattice. One of the reasons is because they are mutually 

] interacting electron systems with both intra- and intersite 
Coulomb repulsions. The former characteristic gives rise to 

] correlation effects and the latter to exchange-correlation ef- 
fects, and both play essential roles in the interesting physics of 

\ these systems. ' In addition, they are very pure systems so their 
observable properties are not caused by impurities, which is 

] an aspect that facilitates the extraction of pure effects that are 
due to mutual interactions between electrons.^- 

\ Among the various phenomena characteristic of organic 
conductors, Wigner-type charge ordering (CO) is particularly 

" interesting because it originates from intersite Coulomb re- 
pulsion.' Surprisingly, recent '■'C-NMR measurements of 0- 

; ET2MM'(SCN)4 (ET = BEDT-TTF, M = Rb, Cs, M' = Zn, 
Co) have revealed the existence of a "glassy" state in quasi- 

] two-dimensional (quasi-2D) organic conductors on the verge 
of CO, which is a spatially inhomogeneous state characterized 

" by a slow charge relaxation with non-Arrhenius behavior.'* X- 
ray diffraction measurements have indicated that this glassy 

" state has an "intermediate-range order", namely, the correla- 
tion length develops well but remains finite in the range of 

[ tens of nanometers,^ which is inseparably connected to giant 
nonlinear conduction, which may act as a thyristor inverter in 
organic electronics.^ 

Because of the pureness of organic conductors, we assume 
that their glassy state does not come from quenched disorder 
but from the mutual interaction between electrons. Although 
it is not well clarified, some eff'orts to identify this glassy state 
have been made on the basis of the 2D extended Hubbard 
model (EHM) with not only intrasite Coulomb repulsion but 
also intersite Coulomb repulsion.^ 

Apart from these eff'orts, there exists a theory of glassi- 
ness that does not rely on the presence of quenched disor- 
der but is based on phase separation (PS) that is frustrated 
by the long-range Coulomb interaction.^'^ Examples of mutu- 
ally interacting electron systems that strongly tend toward PS 
include the t-J model, double-exchange model," and low- 
density electron gases. Theories involving Coulomb frus- 



trated PS'^ have been applied to several classes of materials 
such as high-Tc cuprates,''' colossal magnetoresistive man- 
ganites,'^ and supercritical-fluid alkali metals,"' but have yet 
to be applied to ^-ET organic conductors. 

In this Letter, we report our theoretical results for the 
charge compressibility k. These results are obtained by apply- 
ing the fluctuation-exchange (FLEX) approximation to a sim- 
ple model of 0-ET organic conductors, namely, a 3/4-filled 
2D EHM with both on-site and nearest-neighbor Coulomb 
repulsion. We show that, in a quantum critical regime dom- 
inated by quasiclassical CO fluctuations, /c"' decreases with 
decreasing temperature and eventually becomes negative, 
which indicates that this model has a strong tendency for PS at 
nonzero temperature. By using our numerical results for /c, we 
estimate the correlation length ^ and the period of the elec- 
tron density modulation induced by the long-range Coulomb 
interaction that frustrates the PS. For 0-ET organic conductors 
on the verge of CO, we find that the ratio ^/Zn, can be larger 
than the critical value close to 2, indicating the emergence of 
a glassy state, i.e., a spatially inhomogeneous state charac- 
terized by an extremely slow relaxation and an intermediate 
length scale. 

Our EHM Hamiltonian is defined on a 2D square lattice 
with a lattice constant al as 

i7 = f ^ {c\^Cj^ + H.c.) + f/^ «,|n,x + hihj. (1) 

Here, t is the hopping integral between the nearest-neighbor 
sites, which are denoted by (/, j), c]^ (c„j-) is the creation (an- 
nihilation) operator of an electron on site ; with a spin cr - 
t or i, hicr - c'. Cicr - njl with « being the mean electron 
number per site, «, = n,| -H n,!, and U and V are the on-site 
and nearest-neighbor Coulomb repulsions, respectively. We 
use flL = 1 and t - \ unless otherwise specified. 

We calculate the charge compressibility k = n^^id^ldnY"^ 
for the Hamiltonian eq. (1) at the temperature T . Here, the 
chemical potential yu s ij.{n) is obtained as a function of n 
using 

« = 2 J^G(fc)e'""', (2) 
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where 77 is a positive infinitesimal and ^ denotes 

rXitJi- f^nf^n^^/^^''^^^ ^^'■^ ^' which is a notation that com- 
bines the wave vector k and the fermion Matsubara frequency 
iwi. The single-particle Green's function G(k) is related to the 
self-energy I,{k) through the Dyson equation according to 



G{k) = [iaj,+ti-s,--L(k)]- 



(3) 



where denotes the noninteracting band dispersion. 

For an explicit calculation of k, we adopt the FLEX ap- 
proximation'^' without contributions of the particle-hole 
and particle-particle ladder diagrams, which is the simplest 
conserving approximation for E that includes the nontrivial 
exchange-correlation effect.'^ In this approximation, l,{k) is 
given by 

^(k)^-^ [ G{k-q)Y,Va{q\ (4) 

denotes T J^J^ d9/(27r)- with q, which is a 
notation combining the wave vector q and the boson Mat- 
subara frequency iw^. Here, Ya(<i) represents the exchange- 
correlation interaction potential for a = c, s, as 



= v„{q)l\\ H- v„(^rk(?)], 



(5) 



with vj^q) = U + 4y(cos q^ + cos q^), v^iq) - -U, and;^'(^) = 
- G{k + q)G{k). By combining eqs. (2)-(5), we arrive at a 
self-consistent determination of G for a given n with a relative 
precision of 10"*. Hereinafter, we set U - 6. 

In the FLEX approximation, CO instability is determined 
by the divergence of Vdq, 0). It is thus useful to define the 
"distance" from the CO transition point as 



1 + vc(ecok(eco,o). 



(6) 



where Qco is the wave vector specifying the CO pattern at 
which |yc(?,0)| is maximized. In Fig. 1, we show the con- 
tour curves of Sc on a T-V phase diagram for n - 3/2, where 
the CO transition curve is determined by 6^ - 0. Along 
the transition curve in Fig. 1, we consistently find instabil- 
ity to the unique pattern of checker-board-type CO specified 
by Qco - i^,^) in accordance with previous studies.*^" Be- 
cause Qco is in disagreement with a maximum position of 
Xiq, 0) and is equated to the minimum position of vdq), the 
present CO instability is caused by V, which obliges electrons 
in nearest-neighbor sites to avoid each other in real space. ' 
This disagreement is a generic feature of CO and has been 
discussed for models with several different CO patterns in dif- 
ferent contexts..^' 

From Fig. 1, we see that the intervals between the vari- 
ous contour curves become much larger for a small and 
a nonzero T; typically, Sc < 0.05 and T > 0.2. Such in- 
creasing intervals are an essential feature of critical phenom- 
ena in itinerant electron systems and originate from the self- 
consistent determination of 6c, which includes a feedback ef- 
fect for fluctuations.^^ This self-consistency is partially taken 
into account through the determination of G in the present 
FLEX approximation. Note that the region in which the inter- 
val enlargement occurs roughly corresponds to a quasiclassi- 
cal regime governed by non-Fermi-liquid behavior in the the- 
ory of quantum critical phenomena.^^ In fact, for the EHM, 
the non-Fermi-liquid behavior and its related phenomena oc- 
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Fig. 1. (Color online) A T-V phase diagram for the EHM in the FLEX 
approximation with contour curves for Sc = 0.001, 0.01, 0.05, and 0. 1 (see 
legend for syiTibols) for n = 3/2 and U = 6. 
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Fig. 2. (Color online) Chemical potential Afi = ^(n) - ,u(3/2) as a function 
of natV = 2.2 for T = 0.72 (triangles), T = 0.68 (circles), and T = 0.64 
(squares). The tangent of each curve at n = 3/2 is shown by the solid lines. 



cur in regions dominated by quasiclassical CO fluctuations.^'' 
In Fig. 2, we plot A/i s - i-i(3/2) as a function of n 
for T = 0.72, T = 0.68, and T = 0.64 with V = 2.2. For 
T = 0.72, Afi is a monotonic function with a positive slope 
at « = 3/2. Upon decreasing the temperature to T = 0.68, 
however, the slope at n = 3/2 decreases, and at T = 0.64, it 
is negative. The negative slope of the fi-n curve corresponds 
to a negative k because k = n^^{dfj./dny^ < 0. In Fig. 3, we 
present the entire region of a- < on the same T-V phase 
diagram as shown in Fig. 1, and we find that this region is in 
the quasiclassical regime mentioned in the last paragraph for 
quantum critical CO phenomena. Here, we emphasize that k 
becomes negative for a nonzero T even though it is positive at 
r = 0. 

To clarify the reason for k < 0, we recall that the inter- 
acting Fermi surface is determined by yu = s/c^ + I.{kp, 0) at 
r = 0, where kp is the Fermi wave vector. A similar equa- 
tion is available even for low but finite temperatures; fi is then 
approximated using p. ^ Eko + ReE'^(^o, 0), where ^0 = 9o/2 
with qo being the wave vector at which xiq, 0) is maximized 
and E'*(^o, 0) is the analytic continuation of E(^o, ^oju) with 
(ji)k > onto the origin of the complex frequency plane. There- 
fore, we can attribute the negative a: to a rapid decrease in 
Re2'^(^o, 0) with increasing «, which overcomes the increase 
in 

On the other hand, from eq. (4), we can see that, when 
5c — > and r 0, ReS''(^o> 0) has a positive contribution 
in proportion to ReG'*(^o - 2co> 0) ln(l/5c)- Therefore, we 
can attribute the rapid decrease in ReS''(^^o> 0) to the increase 
in 6c with increasing « in a quantum critical regime. This in- 
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Fig. 3. (Color online) Region of /c < (shaded area indicated by arrow) on 
the same T-V phase diagram as in Fig. 1. 



crease in 6c results from the combined effects of the intersite 
Coulomb repulsion, which makes the position of Qco insen- 
sitive to variations in «, and the existence of the Fermi surface 
with a diameter smaller than l2col, as explained below. 

Because Qco - i^,^) is not sensitive to slight variations 
in n from « = 3/2 for which - 6^ as defined by 

eq. (6) depends on n only through xiQco^^)- Thus, because 
l^ol corresponds to the diameter of the Fermi surface and de- 
creases smaller as n increases, we find that ^o, i-c, the peak 
position ofxiq, 0), moves away from Qco with increasing n. 
This variation in reasonably reduces xiQco, 0), so that 6^ 
increases with n into an off-critical regime, which leads to a 
negative k. Since the discussion to this point does not depend 
on the details of the model, a negative k is expected in general 
for models that display a CO transition specified by a Qco 
whose magnitude is large compared with the diameter of the 
Fermi surface. 

It is important to comment on k for the EHM in the limit 
of y = 0, i.e., the Hubbard model (HM), in which spin- 
density-wave (SDW) fluctuations develop near the half fill- 
ing (n - 1). By using the quantum Monte-Carlo method^^ 
and FLEX approximation,'^ it has been shown that k is al- 
ways positive for the 2D HM even though it tends to diverge 
for n — > 1. In the FLEX approximation, eq. (4) leads to 
the fact that ReS'*(^o, 0) has a contribution proportional to 
ReG'^(A;o - QsDW,0)ln(l/6,) for 6^ ^ with T ^ 0, where 
(5s measures the "distance" from the SDW quantum critical 
point. It is, however, noted that the SDW wave vector Qsdw is 
the nesting vector of the Fermi surface, i.e., gsDW = 2^o and 
then ReG'^(A:o -gsDW, 0) = ReG^{-ko, 0) ~ 0. Therefore, the 
singular contribution oc ln(l/i5s) almost vanishes. This is the 
main reason why SDW fluctuations do not lead to a negative 
K, while CO fluctuations do. 

For the EHM defined by eq. (1), a negative k provides a 
short-range interaction to promote the spatial electron den- 
sity fluctuations 6ii(r), which implies a macroscopic PS due 
to spinodal decomposition. This macroscopic PS is, however, 
cured because the short-range interaction competes with the 
long-range Coulomb interaction between 6n{r) promoted by 
the short-range interaction. This type of arrested PS is called 
Coulomb frustrated PS. 

To see how our results for k within the EHM can lead to 
Coulomb frustrated PS, we first consider the following ex- 
pansion of the free energy with respect to 6n{r): 



SF 



2Df/ 



dr 



CI 



drdr 



, 5n{r)6n{r') 
■-r'\ ' 



(7) 



where Dp and /cp denote the density of states and the 
compressibility for a noninteracting system, respectively. In 
eq. (7), we consider the first and second terms to be the free- 
energy variation for the EHM due to the slowly varying spa- 
tial fluctuations. The last term (with an electronic charge e 
and a dielectric constant e) is taken to be the free energy 
due to long-range Coulomb interaction that frustrates the PS 
and is absent in the EHM. By using the Fourier transform 
n{q) - J (5n(r)e"'' '', eq. (7) can be rewritten as 



SF 



—y 



K 



,|2,2 



n{q)n{-q), (8) 



where the effective Bohr radius is a* 
tity in square brackets in eq. 
mialf^^ at l^rl = K/^)-'/' 



e/(ne^Dp). The quan- 
(8) is minimized for k^/k + 
+ 0. Thus, the system be- 



comes unstable against the spontaneous spatial modulation 
of n{q) + with \q\ = (fl^/;^)"'^^ for k^Ik < -3(/o/fl*)-^^ 
if there is no coupling between fluctuations associated with 
this phase transition. Note that our FLEX calculation of Ik 
for the EHM does not include any effect of these fluctuations 
since they are induced by the competition between the PS and 
long-range Coulomb interaction. 

By adding the quartic term to eq. (7), we can study the 
effects of the coupling between these fluctuations. Defining 
ip{r) = {IqI -\/Dp)6n{r), we then consider the Hamiltonian with 
the coupling constant u as 



4/ 



dr 



Q_ 
4n 



I 



drdr 



■r'\ 



(9) 



where ro 



Kf/K and Q - llil^a*^). In the mean-field 
approximation for eq. (9), the Fourier transform of the cor- 
relation function Q{r - r') = T^^{(p{r)ip(r')} is given by 

- i^r + \q\~ + Q/\q\^ where the parameter r must be 
determined self-consistently by 



d-q 
J2n) 



(10) 



Then Q{q) has a peak at approximately the modulation wave 
number 1^1 - q^= ■\/-r/3 and can be approximated as 



1^1 



qmP 



'-2 



\q\'-3qi\q\ + Q {\q\-qm)^+e 
where the inverse correlation length ^ ' is given by 

r' = ^o' ^2 [doK) - (qM'] /(3qM. 

Substituting eq. (11) into eq. (10), we can obtain 
form 

A T 
2ti t 





3A T 


A/o 


-3{qJof - ^jln 


K 


[(2Zo/fl*)'/3j 



(11) 

(12) 
in the 

(13) 



where A - {jil3>)ul^t and A is the cutoff wave number. Be- 
cause the right-hand sides of eqs. (12) and (13) are equal to 
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Fig. 4. (Color online) f //„ versus T for Vjt = 2.2 and U/t = 6 with t = 
0.108 eV assuming that A = 1 and /o/flg = 1 (flg = 2.56A). The inset 
shows the 2D an'angement of the ET molecules with t, a, and c indicated. 

each other, we have an equation for determining = ■^-r/3, 
and then the period of modulation is given by 1^ = Inlq^- Us- 
ing the solution of this equation for eq. (13), we can obtain the 
correlation length ^. Note that ^ thus obtained never diverges 
unless T -Q; namely, the transition into the modulated phase 
does not occur even though k-pIk < —'iilola*^)^^^ owing to the 
feedback effect for isotropic fluctuations peaked at the modu- 
lation wave number q^^ + 0.^^ 

In this situation, Schmalian et al. have shown that glassi- 
ness emerges instead of the spontaneous density modulation 
if the ratio ^/Zn, is larger than the critical value close to 2.** '^ 
This emergence of a glassy state is signaled by the occurrence 
of an exponentially large number of metastable states and 
long-time correlations characterized by the correlation func- 
tion T{r - r') - r^' \\vs\t^^{ip{r, t)ip{r' , 0)) + 0, leading to an 
extremely slow relaxation and an intermedeate length scale of 
^. Therefore, we can demonstrate that our system has glassy 
dynamics by showing that ^jl^ > 2. 

By using our numerical result of k^/k for the 3/4-filled 
EHM with V/t = 2.2 and U/t = 6, we can actually eval- 
uate ^/Zn, for quasi-2D organic conductors on the verge of 
CO (with the recovery of al and t). For 6'-ET2CsZn(SCN)4, 
the transfer integral and lattice constants shown in the inset 
of Fig. 4 are f = 0.108 eV, a = 9.82 A, and c = 4.87 A.^^ 
We then obtain e = 11 from V = e^/(erNN) = 2.2f and 
Dp = 0.247/(ta^) for the square lattice (^nn - Vo^ + f? and 
flL = ^acl2). The effective Bohr radius is then evaluated as 
fl* = e/ine^Dp) - 2.56A. The cutoff" wave number A should 
be comparable to the magnitude of the reciprocal lattice vec- 
tor and we choose A - 2 ^fln/aL. Numerical values of A and 
lo I flg are expected to be on the order of 1 . Here, we take A - \ 
and /o/flg = 1 as a rough estimation of ^//m- By inserting 
these values together with our results for /cp/zc into eq. (13), 
we obtain ^//m as a function of T, which is shown in Fig. 4. 
From this figure and Fig. 1, we see that ^/Zm has a peak at ap- 
proximately T - 320 [K] at which 6c takes the smallest value. 
We find that ^jl^ is larger than 2, i.e., the value of criterion 
for glassiness, for a broad range of T around this temperature. 
Note that in our calculation, /,„ is on the order of nanometers, 
so that ^ can be over hundreds of nanometers, which seems to 
be too large for the experimentally observed length scale of 
the intermediate-range order 

In actual quasi-2D organic conductors, because of the 
screening by electrons in other layers, the intralayer Coulomb 
interaction is reduced to a screened Coulomb interaction with 



the screening length nearly equal to the interlayer distance 
br^ Therefore, the above values of ^/Zm may be overesti- 
mated. It is, however, noted that the formation of a glassy 
state is virtually unchanged, even considering the screened 
Coulomb interaction instead of the long-range Coulomb in- 
teraction, provided that the condition (Z^ag)'^^ «: 1^ is satis- 
fied.-^ This condition is, in fact, satisfied for the 0-ET organic 
conductors; for example, (/^a*)'^-' ~ a* = 2.56A is much 
smaller than k~ 43.4A-^ for 6i-ET2CsZn(SCN)4. Hence, 
our mechanism can produce a glassy state for actual 0-ET or- 
ganic conductors on the verge of CO. 
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